Domain wall attraction and repulsion during spin-torque-induced coherent motion 
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We calculate the interaction between two magnetic domain walls during their current-induced 
motion. This interaction produces a separation-dependent resistance and also a differential velocity, 
causing domains in motion to experience an effective attraction at large separations and an effective 
repulsion at short separations. In an intermediate range of currents the two domain walls will reach 
a natural equilibrium spacing that depends on the magnitude of the current flowing through the 
material. 



The effects of spin transport in inhomogeneous mag- 
netic systems have important implications for both the 
understanding of fundamental physics and the devel- 
opment of potential applications. Electrical generation 
of spin torque, which is a direct manifestation of the 
conservation of the angular momentum associated with 
spin, permits fast, localized electrical switching of mag- 
netic domains [IHH|, electrical driving of ferromagnetic 
resonance [1HZ], and controlled generation of coherent 
magnons[8J 5j. The study of the motion of domain walls 
induced by this spin torque [TUHH] may lead to novel spin 
torque devices [T9ll21| . The treatment of multiple domain 
walls in a single system[22l [23] is a necessary next step, 
as many of these future devices will require the manipu- 
lation of more than one domain wall at a time. 

Here we calculate spin transport properties and spin 
torque for a pair of n walls in a ferromagnetic semicon- 
ductor separated by a domain of variable size, with its 
magnetization oriented 180 degrees from the magneti- 
zation orientation of the far left and right leads. Us- 
ing a model Hamiltonian[2I] for the coherent transport 
of spin-polarized carriers through a domain wall in the 
absence of spin-orbit interaction, and a piecewise linear 
approximation [25] for the rotation of magnetization in- 
side the domain walls, we calculate transmission and re- 
flection coefficients, and the spin torque as a function of 
separation for each domain wall individually and for the 
system as a whole. We find that the spin torque in each 
domain wall has a distinct non-trivial dependence on the 
separation of the walls. The domain walls repel each 
other when they are very close together and attract each 
other when they are far apart. This suggests that domain 
walls in motion may reach a natural, current-dependent, 
equilibrium distance and has possible implications for the 
motion of multiple domain walls along a "racetrack" [5D] . 

Schematics of the two-wall system are shown in Fig. [T] 
Two Neel type domain walls are separated by a domain, 
the magnetization of which is antiparallel to the left and 
right leads. The domain walls are oriented such that 
they represent a full 360 degree rotation in magnetiza- 
tion, rather than 180 degrees and back. Spin polarized 
carriers are injected from the left lead. The exchange 
field for such a system can be approximated as: 



B = £> [sin Q(x)$l + cos 9(x)z] , 



(1) 



where the form of 9{x) is shown in Fig. [2] for a pair of 
1.8 nm domain walls separated by 5 nm, and is based 
on a realistic form for the magnetization inside a domain 
wallHaiMl. 
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FIG. 1. (color online) Schematic representation of two Neel 
type 7r domain walls separated by a domain with magnetiza- 
tion antiparallel to the magnetization in the leads. 



The Schodinger equation for the system with this ex- 
change field is 
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where A is the energy splitting between carriers of oppo- 
site spin orientation in the ferromagnetic material. 

This Schrodinger equation cannot be analytically 
solved for the 9{x) shown in Fig. [2] We approximate 
the magnetization inside the domain walls as a piecewise 
linear function, such that 6i{x) — (4>i(x — Xi))/di, where 
4>i is the total magnetization rotation in segment i, Xi is 
the leftmost position of the segment, and di is the width 
of the segment. For each linear domain wall segment, we 
can solve Eq. [2] by transforming to a rotating basis|27). 
The rotation matrix: 



Ri = e 



cos 



sin 



(3) 



defines tpi — Ri^Pi and removes the 9 dependence from 
the off-diagonal potential matrix: 



R7 1 
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- cos 6i 



Ri = cr. 



(4) 
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and yields a modified Schrodinger equation: 
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This Schodinger equation can be solved analytically for 
the wavefunctions in each piecewise-linear segment of the 
domain wall. We use transfer matrices to connect the 
wavefunctions in each linear segment, and then solve for 
the total reflection and transmission coefficients, with 
and without spin flip, as well as the wavefunctions for 
the central domain. 
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FIG. 2. (color online) Magnetization as a function of position 
for a two-domain wall system. The domain walls here are 1.8 
nm and the separation between them is 5 nm. The blue boxes 
show the position of the domain walls. 

After obtaining the full wavefunctions for the entire 
system, we define a spin current density[28 : 
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2im 
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The tensor Q has a flow direction in real space as well 
as a direction in spin space. As our transport model is 
one-dimensional, the real-space flow direction lies solely 
along the x direction, and we write Q as a vector with 
components corresponding to the appropriate spin-space 
directions. As this spin current is not a conserved quan- 
tity, we can then define the spin torque per unit area as 
the amount of spin current lost to the domain wall during 
transport [2*5]: 



N 



DW 



Qt-Qi 



(7) 



The total spin torque is then calculated by integrating 
the transmission and reflection coefficients over the car- 
rier population. Calculations shown here are for a pair of 
1.8 nm domain walls in a material with a spin splitting of 
100 meV and an effective carrier mass of 0.45 m e where 
m e is the mass of the bare electron. The calculations are 
performed with a temperature of 110 K, and the carrier 
density is calculated to be ~ 10 19 cm~ 3 . These param- 
eters are representative of GaMnAs (without spin-orbit 
interaction), but the general qualitative features found 
here do not depend on the detailed quantitative param- 
eters chosen. 



Fig. [3] shows calculated probabilities for transmission 
and reflection of the carriers, with and without spin flip, 
for these two domain walls, for four domain wall separa- 
tions. At a distance of 10 nm (a), we see a transmission 
spectrum with a number of resonance peaks in the energy 
region below the spin- splitting A, and largely transmis- 
sion without any spin flip above A. This behavior is 
that of a spin-dcpcndcnt double barrier resonant tunnel- 
ing structure. At a distance of 5 nm (b), the resonant 
behavior remains, with fewer peaks. At a distance of 1 
nm (c) , the peaks are even fewer, and there is additional 
transmission below A as the domain walls are getting 
close enough to interact more with each other. When 
the two walls are brought into contact with each other 
(d), the transmission spectrum has many of the features 
of a single 2ir wall (inset [25 ). 
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FIG. 3. (color online) Probabilities for transmission and re- 
flection with and without spin flip, for a pair of 1.8 nm Neel 7r 
walls with separation between them of (a) 10 nm, (b) 5 nm, (c) 
1 nm, and (d) nm. Inset: The same transmission/reflection 
spectrum for a 3.6 nm 2n wall. 

The calculated energy-resolved components (torkance) 
of the total spin torque acting on both walls from Eq. (*7]) 
are shown in Fig. [4] We identify the spin torque as adi- 
abatic (proportional to VM(r), and thus parallel to x) 
or non-adiabatic (proportional to M(r) x VM(r), par- 
allel to y). We see the same resonant behavior in the 
spin torque that we saw in the transmission spectrum for 
large separations of the domain walls, increasing interac- 
tion with decreasing separation, and spin torque similar 
to that of a 2ir wall when the domain walls are brought 
into contact. 

Fig. [5] shows the total spin torque for the double wall 
system as a whole, as a function of the separation of the 
two walls, and compares it to the spin torque in individ- 
ual 7r and 27T walls. The magnitude of the spin torque 
increases as the two walls move a small distance away 
from each other, and then falls to a saturation level that 
is higher than the spin torque in a 2-7r wall the same size 
as the two walls together, and higher than the twice the 
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FIG. 4. (color online) Spin torque as a function of energy, for 
a pair of 1.8 nm Neel n walls with separation between them 
of (a) 10 nm, (b) 5 nm, (c) 1 nm, and (d) nm. Inset: The 
same calculation for a 3.6 nm 2-k wall. 
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FIG. 5. (color online) Spin torque through the entire two wall 
system (solid lines). Dashed lines are the same calculated 
value for a single 2ir wall of 3.6 nm. Dotted lines are twice 
the calculated value for a single -k wall of 1.8 nm. 

The spin current from Eq. © is a position-dependent 
quantity, so the torque acting on any region of space, 
produced by carrier motion through the region of space, 
can be determined from Eq. Q. These expressions, and 
the values of the spin-dependent wave functions, have 
been used to obtain the torque acting separately on the 
upstream and the downstream domain walls (relative to 
the current direction). Fig. [6] shows the breakdown of the 
spin torque for each individual domain wall as a function 



of the separation of the two walls. The velocity of one of 
these domain wall can be estimated as 



9fJ± 
M, 



N 



(8) 



where N is the total spin torque in h flips per second per 
cm 2 and M s is the saturation magnetization of the ma- 
terial. As this is directly proportional to the spin torque, 
this picture of the torque represents the motion of each 
domain wall. 
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FIG. 6. (color online) Spin torque (a) along the direction of 
motion (x) and (b) perpendicular to the direction of motion 
(y) as a function of domain wall separation, for an upstream 
(solid line) and a downstream (dashed line) 1.8 nm Neel n 
wall. 

If the motion of the domain wall can be constrained 
so that there isn't any motion in the y direction (such 
as is possible in certain thin films and nanowires), the 
motion of the domain walls is in the x direction, which 
is along the same direction as the current flow. The spin 
torque in Fig. [6]Ja) shows three distinct regimes of inter- 
action between the velocities of the two walls along the 
x direction. 

For very small separations between the two walls, the 
upstream wall moves in the —x direction and the down- 
stream wall is moves in the +x direction, so the separa- 
tion between the walls grows. From the perspective of a 
coordinate system whose origin lies at the midpoint of the 
two walls the two walls appear to be repelling each other. 
As the separation between the walls increases, the veloc- 
ity of the upstream wall slows, stops, and begins to move 
in the +x direction. In this second regime of interaction, 
the two walls move together in the +x direction. As the 
velocity of the downstream wall decreases, this results in 
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a separation where both walls experience equal amounts 
of torque, and thus move with equal velocities. For still 
larger separations the downstream wall eventually slows, 
stops, and begins to move in the —x direction, while the 
upstream wall continues in the +x direction. From the 
perspective of a coordinate system whose origin lies at 
the midpoint of the two walls the two walls now appear 
to be attracting each other. 

The domain walls used for this calculation were 1.8 nm 
wide, as previous calculations showed a high spin torque 
for a 7r wall of this size, with a maximum around 2-3 
nm width [25]. The calculations were repeated for other 
widths of domain walls. For thinner domain walls the 
same patterns of interaction between the two walls are 
found, although the equilibrium separation between up- 
stream and downstream walls was closer together. For 
much thicker domain walls (thicker than around 4-5 nm) , 
the crossover of the velocity of each wall from positive to 
negative as a function of separation no longer occurs, and 
the interaction is repulsive for the entire range of wall sep- 
arations. The upstream wall's velocity approaches zero 
as the separation increases, and the downstream wall's 
velocity saturates to a positive value as the separation 
increases. 

We have shown that the velocities of two separated 7r 
domain walls depend in a sensitive and nonlinear fashion 
on the wall separation. The energy-dependent electronic 
transmission and reflection, and the energy-resolved spin 
torque (torkance) change from those of a spin-dependent 
double barrier resonant tunneling structure (when well 
separated) to those of a 2n wall (when close together). 
Three interaction regimes have been found for the depen- 
dence of the differential spin torque of the two domain 
walls on wall separation — repulsion, motion together, 
and attraction. The existence of these three regimes in- 
dicates a stable equilibrium separation between two do- 
main walls in the presence of current flow; thus two or 
more domain walls could be configured to move with a 
stable separation along a racetrack [2D] with the applica- 
tion of an appropriate current. 
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